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We study the SU(2) gauge-Higgs system in D — 4 dimensions, and analyze the influence of the 
fundamental-representation Higgs field on the vortex content of the gauge field. It is shown that 
center projected Polyakov lines, at low temperature, are finite in the infinite volume limit, which 
means that the center vortex distribution is consistent with color screening. In addition we confirm 
and further investigate the presence of a "Kertesz line" in the strong-coupling region of the phase 
diagram, which we relate to the percolation properties of center vortices. It is shown that this 
Kertesz line separates the gauge-Higgs phase diagram into two regions: a confinement-like region, 
in which center vortices percolate, and a Higgs region, in which they do not. The free energy of the 
gauge-Higgs system, however, is analytic across the Kertesz line. 

PACS numbers: 11.15.Ha,12.38.Aw 



I. INTRODUCTION 

Lattice studies strongly support the vortex picture of 
confinement and the importance of the center degrees of 
freedom (see review Q and references therein). Until 
recently, these numerical investigations concerned pure 
gluonic QCD only. If dynamical fermions are included, 
it is well known that the long range part of the potential 
between static charges changes qualitatively: the string 
breaks and the potential levels off at a screening distance 
rrj. String breaking has been observed numerically not 
only for gauge fields with dynamical fermions, but also 
for gauge-Higgs theories |2| ■ In this work we investigate 
the influence of dynamical matter fields on the distri- 
bution of center vortices. Our model is lattice SU(2) 
gauge-Higgs theory with Higgs fields in the fundamental 
representation. 

The SU(2) gauge-Higgs model is defined by the action 
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where Sw is the usual Wilson plaquette action, and $ = 
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FIG. 1: Schematic phase diagram for the SU(2)-Higgs sys- 
tem at zero temperature (top), and on lattices with a fixed 
extension in the time direction (bottom). 



(0*), with 4>i,4>2 € C, is a massive scalar field in the 
fundamental representation. 

A schematic phase diagram for the SU(2)-Higgs model 
is depicted in Fig.^ ln the "confined" phase, the poten- 
tial between fundamental charges rises linearly at inter- 
mediate distances. Due to color screening of fundamental 
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charges, there is string breaking at some finite distance 
r , and the potential levels off. In the Higgs phase, the 
Higgs mechanism is at work, and the potential is Yukawa- 
like; the string tension vanishes at all separations. How- 
ever, these are not thermodynamically distinct phases. 
The phase diagram is connected, in the sense that one 
can always find a path between any two points in the 
phase diagram which avoids any non-analyticity in ther- 
modynamic quantities. The transition line which might 
have separated the Higgs and confinement phases ends at 
a critical point in the strong coupling region, and from 
there changes over to crossover behaviour. Thus if we 
speak of the "Higgs" or the "confinement" phases in this 
model, we are aware that we speak rather of regions be- 
longing to a single phase of the system Q. 

At finite temperature and K — there is an additional 
phase, namely, the deconfined phase. In this phase the 
quark-gluon plasma arises and fundamental charges are 
set free. For k > there is no true phase transition be- 
tween the "confined" and "deconfined" phases, but only 
crossover behavior indicated by the dashed line in Fig.^ 
(lower figure). Hence we still have only one phase in the 
[3 — k plane. 

The vortex theory of confinement was put forward at 
the end of the 1970s (see references in review [l|); early 
applications to SU(2) gauge-Higgs theory can be found 
in refs. 0, 0- According to the vortex picture, tube- 
like (D=3) or surface-like (D=4) objects, carrying quan- 
tized amounts of magnetic flux, play a crucial role. The 
standard procedure to identify these objects (the center 
vortices) on the lattice is to first gauge-fix the lattice con- 
figurations according to a gauge-fixing condition on link 
variables in the adjoint representation. Then the cen- 
ter degrees of freedom are singled out by projecting each 
SU(2) group-valued link variable to the closest Z 2 center 
element. In this work we use the direct maximal center 
gauge (DMCG) |(| fixed by the over-relaxation method. 
Maximizing (|Tr[C/ M (a;)]| 2 ), DMCG shifts the link vari- 
ables U fl (x) as close as possible towards the center ele- 
ments of SU(2). Center projection in DMCG has proven 
to be a useful tool for isolating the relevant degrees of 
freedom on the lattice, although other adjoint gauges, 
and more sophisticated gauge-fixing techniques, arc also 
in use. Via center projection, the SU(2) lattice is mapped 
to a Z2 gauge field configuration. The excitations of Z 2 
lattice configurations are the P-vortices. P-vortices form 
closed surfaces (in four dimensions) or closed loops (in 
three dimensions) on the dual lattice, and are composed 
of P-plaquettes. A plaquette on the projected lattice is 
said to be pierced by a P-vortex, and is called a "P- 
plaquette" , if the product of the four projected links of 
the plaquette yields the nontrivial element of Z 2 , i.e. —1. 
Center projection in an adjoint gauge can be viewed as 
a tool to locate the position of thick SU(2) vortices from 
the location of P-vortices on the projected lattice. 

Here we are interested in the dependence of the cen- 
ter vortex distribution on couplings and temperature. 
The vortex model at finite temperature, and the be- 



haviour of P-vortices at the transition to the high tem- 
perature phase, have been studied in recent years in 
refs. 0, la El EH]- In those articles it was shown 
that center dominance — i.e. the correspondence of pro- 
jected and unprojected observables sensitive to infrared 
physics — holds also at finite temperatures. It was found 
H> Ei EH E3 that the vortex density p — the frac- 
tion of plaquettes on the projected lattice which are P- 
plaquettes — drops substantially across the deconfinc- 
ment transition and, more importantly, that vortices do 
not percolate in any space-slice, i.e. there is no perco- 
lation of vortex lines in a three-volume with either the 
x, y, or z coordinate held fixed. The absence of vortex 
percolation in a space-slice implies that center projected 
Polyakov lines have non-zero expectation values in the 
deconfined phase l9lllll. 

Recent findings |12l Ha | indicate that the vortex den- 
sity also decreases sharply in the gauge-Higgs model, at 
the transition from the "confined" phase to the Higgs 
phase. This sharp decrease in vortex density occurs even 
in the crossover region at small values of /3, where the 
phase transition line has ended. It was suggested by 
Langfeld |l2j that the cross-over line at small /3, where 
there is a sudden drop in the vortex density, could be 
a "Kertesz line" of the sort found in the Ising model 
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A Kertesz line is a line of percolation transitions 
which is not associated with a thermodynamic transition, 
e.g. from an ordered to a disordered state. Our new cal- 
culations provide further support for the existence of a 
Kertesz line in the gauge-Higgs system. Although there 
is no thermodynamic transition from a "confined" to a 
Higgs phase, and indeed (in contrast to pure gauge the- 
ory) there is no true confined phase in this system, we 
nonetheless find a line of center vortex depercolation. 

II. FINITE TEMPERATURE 

It is most efficient to carry out Monte Carlo simulations 
of the gauge-Higgs system in unitary gauge, where $ = 
(q) , (f> G [0, 00] , and only one degree of freedom has to be 
simulated for the Higgs field. In the unitary gauge we 
have 

Z = J -D[U]V[4>] exp (-SW - Sau) (3) 

X ( 4 ) 
-k2(^(i)0(x+A)'&[C^(*)]) 

f_L . X 

with cj) — V <I>f <£> and integration measure 

[m=l[ [°° d<p(x)<p 3 ( X ). (5) 

J x J 

We used the Metropolis algorithm for the numerical sim- 
ulation of the system on D = 4 dimensional Euclidean 
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lattices of dimensions N t *N^. Data was taken on lattices 
separated by 50 Monte Carlo steps. For the measurement 
of center projected observables, we fixed to DMCG us- 
ing the over relaxation method and then carried out the 
projection. 

At finite temperature, confinement and deconfinement 
can be studied using the Polyakov loop observable 
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and F is the free energy of a single static quark relative 
to vacuum at temperature T. Without dynamical matter 
fields, we have 



lim (\L\) = 



(8) 



in the confined low temperature phase. If matter is in- 
cluded, screening becomes possible. In that case (\L\) is 
non-zero, and F is finite, in the large volume limit. 

Since the idea of the vortex model is to describe the 
confining properties of gauge fields by center degrees 
of freedom, the question arises whether this observed 
screening can also be seen after singling out these vari- 
ables. We have carried out center projection in DMCG 
without (k = 0) and with (k = 0.25) scalar matter fields 
on lattices of size N t * with N t — 4 and N s vary- 
ing from 8 to 28. The quartic coupling has been set to 
A = 0.5, and data has been taken on 1000 configurations 
for each parameter set. Some results are shown in Fig. [21 
for ft = 2.1 and 2.2, where the system is in the low tem- 
perature "confined" phase. 

At the chosen j3 v alue s and k — 0, the expectation 
value of vs. 1/ y/Nf on unprojected lattices extrap- 
olates linearly to zero in the N s — > oo limit of infinite spa- 
tial volume. In contrast, at finite coupling k = 0.25, (\L\) 
tends to a non-zero constant at infinite spatial volume. 
This well known feature, associated with the screen- 
ing of static charges by dynamical matter fields, is also 
found on center projected lattices. Due to suppressed 
self energy contributions and ultraviolet fluctuations in 
the projected configurations, the expectation value of the 
Polyakov loop observable (|£p|) on projected lattices is 
larger than the corresponding expectation value on 
unprojected lattices, but the qualitative behaviour of the 
two observables, in the large- volume limit, is the same. In 
fact, the ratio of the projected and unprojected Polyakov 
expectation values is virtually constant, and nearly inde- 
pendent of both lattice extension N s and coupling k for 
all values of (3 = 2.1, 2.2 and 2.25, as seen in Fig. [3J 

At k = the Polyakov expectation value (\L\) is zero 
in the infinite volume limit at low temperature, and non- 
zero beyond the deconfinement transition. This transi- 
tion can be understood in the framework of the vortex 
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FIG. 2: Expectation values of the Polyakov loop in the "con- 
fined" phase. Open symbols denote measurements at n — 0.0, 
filled symbols at n = 0.25. Squares are used for measurements 
on unprojected, full SU(2) configurations and circles for cen- 
ter projected fields. 
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FIG. 3: Ratio of projected to unprojected Polyakov loop ex- 
pectation values in the "confined" phase. Open symbols de- 
note measurements at k — 0.0, filled symbols at k = 0.25. 
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FIG. 4: Expectation values of the Polyakov loop in the decon- 
fined phase. Open symbols denote measurements at k = 0.0, 
filled symbols at k — 0.25. Squares are used for measure- 
ments on unprojected, full SU(2) configurations and circles 
for center projected fields. 



model HHOilGil. The projected Polyakov line (\L P \) 
has the same transition from zero to non-zero values at 
the deconfinement temperature, as already reported in 
ref. 0. 

The inclusion of a massive scalar field does not lead 
to a qualitative change of the system in the high tem- 
perature phase. This is reflected by our results for (\L\) 
and (\L P \) at (3 = 2.3 and 2.4, seen in Fig.Hand Fig.0 
Again the projected and the unprojected measurements 
are consistent. As in the "confined" phase, the projected 
values are higher than the unprojected ones, but their ra- 
tio is almost constant with respect to N s . There is some k 
dependence in projected and unprojected Polyakov lines 
right above the high T phase transition (which is just 
below (3 = 2.3 at N t — 4), but the ratios of projected 
and unprojected lines are nevertheless almost constant, 
and the k dependence in (|L|) and (\Lp\) is rather small 
at still higher temperatures in the deconfined regime. 

From these calculations we can conclude that we see 
center dominance for Polyakov loops also in the presence 
of scalar matter fields, both above and below the finite 
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FIG. 5: Ratio of projected to unprojected Polyakov loop ex- 
pectation values in the deconfined phase. Open symbols de- 
note measurements at k — 0.0, filled symbols at k = 0.25. 



temperature phase transition. This means that and 
(|Lp|) are in the infinite volume limit cither both zero, or 
both non-zero, at any temperature with or without dy- 
namical matter fields. Once again, the long range physics 
of the system seems to be encoded in the center degrees of 
freedom. In particular, the presence of a dynamical mat- 
ter field must alter the spatial distribution of P-vortices, 
in such a way that P-vortex fluctuations no longer bring 
(|£.p |) to zero in the large volume limit, even in the low 
temperature regime. 



III. ZERO TEMPERATURE 

Our results for T = have been obtained on hypercu- 
bic 10 4 lattices for (3 < 1.9, on 12 4 lattices for ft = 2.1, 
and on 16 4 lattices at j3 = 2.3. For (3 — 0.25, we have 
performed a detailed scan over k on a 16 4 lattice. Each 
measurement has been done at A = 1.0 using 800 config- 
urations. Otherwise, we have used the same parameters 
and methods as for the simulations at finite temperature. 

Due to the small overlap of flux-tube states with 
gluelump states, it is difficult to measure string breaking 
in the "confined" phase using only Wilson loops, since 
the breaking (and hence perimeter-law behavior) occurs 
only for rather large loops, where extremely high statis- 
tics are needed. 1 This problem is reflected in the vortex 
model. Two of us (R.B. and M.F.) have carefully inves- 
tigated the Wilson loop observable in a wide range of the 
phase diagram (F0 < [3 < 2.3,0.0 < n < 1.1) on lat- 
tices up to 22 4 13]. As in the unprojected, full system, 



1 Nevertheless, we note that recently adjoint string breaking has 
been successfully measured using Wilson loops in ref. 111! , using 
a pow erfu l algorithm for noise reduction devised by Liischer and 
Weisz IT7I. 
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FIG. 6: P-vortex density p in the j3 — « plane. Lattice sizes 
are indicated in the text. The phase transition line is taken 
from 0. 



the influence of the Higgs on the Wilson loop observable 
W(R, T) can be effectively described — for our available 
small loops i?, T < 4 — by some shift (3 — ► /3 e //. Center 
dominance holds again. 

At the transition to the Higgs phase, the string ten- 
sion measured by Wilson loops disappears even for small 
distances. We have found that the phase transition line 
in the gauge-Higgs model is closely correlated, even for 
small lattices, with a rapid drop in the vortex density 
p, as shown in Fig. In this figure we plot the den- 
sity of P-vortex plaquettes in percent as contours in the 
(3 — k plane. Data for the phase transition line are taken 
from 0| . In the "confined" phase, the vortex density 
is almost independent of k, whereas the vortex density 
decreases rapidly, with increasing k, in the neighborhood 
of the phase transition line. This rather sudden decrease 
in vortex density with k, at fixed (3 was first observed by 
Langfeld [3, and by two of us (R.B. and M.F.) [l]|. It 
is interesting that this region of rapid decrease in density, 
seen in Fig. El appears to extend beyond the actual ther- 
modynamic transition line, all the way to the (3 — axis. 
We may ask whether this is simply crossover behavior, 
or an indication of some more interesting phenomenon. 



A. The Kertesz Line 

A "Kertesz line", originally discovered in Ising spin 
systems, is a line of percolation transitions which is free 
of any non-analyticity in the corresponding free energy. 
This is therefore not a line of thermodynamic phase tran- 
sitions, as usually defined. It was originally suggested by 
Langfeld that the rapid decrease of P-vortex den- 
sity in the gauge-Higgs theory, in a region where there is 
no thermodynamic transition, might also be associated 
with a Kertesz line of P-vortex percolation transitions. 
We will now report on some results which support that 



FIG. 7: P-vortex density p, weighted average vortex size 
s w , and the one- link gauge-Higgs energy density Ogh for 
j3 — 0.25. Open symbols denote measurements after center 
projection, for filled symbols (labelled with "sm" ) elementary 
vortex cubes are removed in addition. 

idea. 

Our simulations are performed on 16 4 lattices at (3 — 
0.25, which is far below the end of the thermodynamic 
phase transition line around (3 = 1. In Fig. [7] we display 
the following observables: 

• The one-link contribution to the gauge-Higgs en- 
ergy density 

O gh := (Re [<S>\x)U^x)<P(x + A)] >. (9) 
This is a local, thermodynamic observable. 

• The P-vortex density p. 

• The weighted average vortex size s w . This is the 
relative size ||ci( )||/JVp of the vortex Cj( a ) contain- 
ing the P-plaquette P a , averaged over all Np P- 
plaquettes in the projected gauge field configura- 
tion 2 : 




Here ||cj|| is the number of plaquettes in vortex 
Cj, i(a) the index of the vortex containing plaque- 
tte P a , and N v is the number of vortices in the 
projected configuration. In the limit of infinite lat- 
tice volume and finite vortex density, s w = 1 if all 



2 This observable resembles the average cluster size S of percola- 
tion theory. In our case the clusters are the vortices. S differs 
from our s w in that i) it is not the average absolute size, but the 
average fraction with respect to all P-plaquettes; and ii) all clus- 
ters resp. vortices, including the percolating one, are included in 
the sum. 
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P-plaquettes belong to a single percolating vortex, 
and s w = for the depercolation case. This is a 
highly non-local observable, and it is the quantity 
which we use to detect a percolation transition. 

At j3 = 0.25, it can be seen from Fig. that the one- 
link energy density Ogh is a smooth function of k, and 
displays no evidence of a transition. The vortex density, 
although dropping rapidly in the interval k € [0.6,0.9], 
also shows no sign of any discontinuity. 

The vortex observables in Fig. have been computed 
in two ways. Open symbols denote measurements on cen- 
ter projected configurations. For filled symbols (labelled 
with "sm"), the first s tep of a vortex smoothing proce- 
dure, described in ref. [lj, has been applied. This first 
step completely removes the smallest (one-cube) vortices 
consisting of six P-plaquettes. We have checked that 
down to (3 = 1.0 this removal does not alter projected 
Creutz ratios x(R,R), R > 2. The smoothing procedure 
decreases the vortex density p depicted in Fig. [7| slightly, 
but does not qualitatively change the shape of the curve. 
Smoothed or not, the vortex density decreases continu- 
ously from just under 50% to zero as k increases, with 
no evidence of any sudden change in phase. 

There is, however, a clear sign of a sudden transition, 
around k = 0.74, in the weighted average vortex size 
s w . The plot of s w in Fig. shows a sharp deperco- 
lation transition for this observable. For n < 0.7, the 
majority of P-plaquettes belong to one big, percolating 
vortex surface, but for k > 0.74, this surface is split into 
many small vortices. This is even more pronounced if 
we remove the smallest (one cube) vortices using the 
smoothing procedure. After smoothing, but below the 
transition, the average vortex contains over 90% of all P- 
plaquettes. Just after the transition, this number drops 
to under 1.5%. This is a fairly convincing sign of the ex- 
istence of a percolation transition. We have also checked 
for percolation directly, by measuring the spatial exten- 
sion of the largest vortex surface on the lattice. In the 
Higgs phase, the largest vortex always fits inside a hy- 
percube which is smaller than the full lattice, while in 
the "confined" phase this is not possible; the largest con- 
nected vortex surface extends through the entire lattice, 
irrespective of lattice size. Thus we find center vortex 
percolation in the "confined" phase of the gauge-Higgs 
theory, and no percolation in the Higgs phase. There is 
a first-order phase transition line, which has an endpoint 
in the interior of the (3 — k plane but then continues as a 
Kertesz line, completely separating these two regions of 
the phase diagram. 

It is quite natural that there should be no vortex per- 
colation in the Higgs phase, since percolation is a neces- 
sary condition for achieving an area law for Wilson loops, 
which of course is absent in this phase. But an asymp- 
totic string tension is also absent in the "confined" phase, 
due to color screening, yet in this case we do have vortex 
percolation. This is not a paradox, however. The vortex 
confinement mechanism requires that vortex piercings of 
the minimal area of a large Wilson loop are uncorrelated, 



and percolation alone is not enough to ensure this prop- 
erty. As an example, consider a percolating P-vortex 
surface having the form of a branched polymer (a form 
which often arises in numerical simulations of random 
surfaces). In this case each piercing of a plane surface 
will be accompanied by a second nearby piercing, with 
the pair contributing no net center flux, and therefore 
no net disordering, to the Wilson loop. In this example 
there is no confinement, even though the P-vortex per- 
colates throughout the lattice. Percolation is therefore a 
necessary, but not a sufficient condition for confinement 

G3- 



IV. CONCLUSIONS 

We have observed center dominance in SU(2) gauge- 
Higgs theory for the Polyakov line observable, both at low 
and high temperatures, throughout the f3 — k coupling 
plane. Polyakov line expectation values and (|ip|), 
on the full and projected lattices respectively, arc cither 
both zero, or both non-zero, in the infinite volume limit. 

There is no true confining phase in a gauge-Higgs the- 
ory at k > 0, or in any gauge theory with matter in the 
fundamental representation of the gauge group. Never- 
theless, we have found a non-local observable, sensitive 
to vortex percolation, which distinguishes between the 
confinement-like and Higgs regions of the phase diagram. 
The confinement-like region is characterized by an area- 
law decay of Wilson loops up to a certain string-breaking 
scale, while in the Higgs-like region there is no area-law 
falloff at any scale. The two regions are separated by 
a line of thermodynamic first order transitions, which 
turns into a Kertesz line at j3 < 1. As in the Ising spin 
system, the Kertesz line is a line of percolation transi- 
tions, with the free energy analytic across the transition. 
In the confinement-like region, center vortices percolate 
throughout the lattice. In the Higgs region, they do not. 

The sharp transition between the Higgs and 
confinement-like regions has been seen in other ways. 
Langfeld has noted that after fixing to Landau gauge, 
there is remnant unfixed global symmetry, and that this 
symmetry is unbroken in the confinement-like phase, and 
broken across the Kertesz line in the Higgs phase. A sim- 
ilar observation, this time in Coulomb gauge, has been 
made very recently b y tw o of us (J.G. and S.O.) and 
Zwanziger in ref. |2(J, |21| . It should be noted that the 
order parameter for remnant symmetry breaking, if ex- 
pressed as a gauge-invariant observable, is highly non- 
local, as is the order parameter for percolation. 

The fact that vortices percolate in the confinement-like 
region of the gauge-Higgs phase diagram, yet the asymp- 
totic string tension is zero, demonstrates that while vor- 
tex percolation is a necessary condition for confinement, 
it is not also a sufficient condition. Vortex percolation 
without confinement implies that vortex piercings of a 
plane are not entirely random, but are paired in some 
way. An example of such a percolating but non-confining 
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surface is a branched polymer, but the actual structure of 
percolating center vortices, in the confinement-like region 
of gauge-Higgs theory, is not yet known. 
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